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Abstract. Recent ROSAT studies of narrow-line Seyfert 
1 galaxies revealed these objects (with FWHM Hp < 
2000 km s~ 1 ) generally show steeper soft X-ray spectra 
than broad line Seyfert l's and that there are no AGN with 
broad lines and steep X-ray spectra. We derive a simple 
theoretical model which explains this observed correlation 
between the line width and the spectral index for Seyfert 
1 galaxies. Assuming the line width is due to gravitational 
velocity dispersion, it is determined by the radius of the 
broad line region. Sources with steep X-ray spectra (for 
a given luminosity) have a stronger ionizing power than 
flat-spectrum sources with comparable luminosity, which 
implies that the BLR is formed at relatively larger dis- 
tances from the central source, and hence has a smaller 
velocity dispersion and a smaller observed FWHM. We 
test the model over a hetrogeneous (normal and narrow- 
line) sample of some 50 AGN finding a good agreement 
with the data. 

Key words: Galaxies: nuclei — Galaxies: active — X- 
rays: galaxies — Black hole physics 



1. Introduction 

Boiler, Brand & Fink 1996 (hereafter BBF96) report the 
observation with ROSAT of a sample of 49 Narrow-Line 
Seyfert 1 galaxies (NLS1), finding that NLS1 have gener- 
ally steeper soft X-ray continua than normal Seyfert Is. 
(However, not all NLS1 are remarkably steep). As a re- 
sult, looking in the FWHM-a plane, there are no AGN 
with broad lines and steep X-ray spectra. 

BBF96 tested NLS1 models assuming they are Seyfert 
1 galaxies with pole-on orientation, warm absorption, 
thick BLR and smaller black hole mass and/or higher ac- 
cretion rates. The latter showed some promise, but have 
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drawbacks as well, as some NLSls do not have steep spec- 
tra, and on the other hand a smaller black hole mass does 
not necessarily produce narrower lines (see Sect. 2). 

We suggest a simple physical explanation to the ob- 
served correlation between the soft X-ray spectral index 
and the H/3 line width. Assuming the narrower permitted 
lines reflects a lower Keplerian velocity, we explain the 
observed correlation by showing that a steeper spectrum 
has a stronger ionizing power, and hence the BLR form 
at a relatively larger distance from the central source. (In- 
deed, the characteristic distance of the BLR estimated by 
this method, agrees well with the BLR sizes determined 
by reverberation mapping (Wandel 1996;1997)). 

Below we formulate this scanario by combining the as- 
sumptions of (1) Keplerian velocity Doppler width for the 
broad lines, (2) a power-law spectrum of the ionizing con- 
tinuum luminosity, and (3) characteristic values for the 
ionization parameter and physical conditions in the BLR. 



2. The analytic approach 

We derive a simple analytic relation between the width 
of broad emission lines and the continuum slope of the 
central engine. The observed velocity dispersion is deter- 
mined by the mass of the black hole and the Broad line 
region radius. (Note however that the BH mass does not 
solely determine the BLR velocity; for example, in Bondi 
accretion we have M ~ R 2 where R a ~ GM/kT is the 
accretion radius, so that M ~ M 2 . If the efficiency is con- 
stant, i ~ M 2 than R ~ is ~ M so v oc M/R ~ const). 
Our basic idea is as follows: Sources with steep X-ray spec- 
tra (for a given luminosity) have a stronger ionizing power 
than sources with flat spectra. Since the ionization param- 
eter at the BLR has a characteristic value, a larger ionizing 
power implies that the BLR is formed at a relatively larger 
distance from the central source, and hence has a smaller 
velocity dispersion and a smaller observed FWHM. 

We make the following assumptions (rather standard 
in the literature) : 
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Fig. 1. X-ray continuum slope, obtained from a power-law fit, 
versus FWHM of H/3 for Seyfert 1 galaxies (taken from Walter 
& Fink 1993) and NLS1 (BBF96). Lines of constant X-ray 
luminosity are overplotted for the relation given by Eq. 7. 
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a. The width of the broad lines is induced by Keplerian 
motion in the gravitational potential of the central mass. 
The full width at half maximum is given by: 



FWHM 



(i) 



where M is the mass of the central black hole and R the 
radius of the broad line region. 

b. The physical conditions in the ionized gas emitting 
the broad lines are characterized by the ionization parame- 
ter U, the ratio of ionizing photons to electrons (cf. Netzer 
1990) defined by 

U = Q ton /4:irR 2 cn e where Q ion = J^ g l(E)^§- is the 
ionizing photon flux (number of ionizing photons per unit 
time) 1(E) is the monochromatic luminosity of the central 
source, per unit energy, and n e is the electron density. 

The radius of the BLR may be written as 



R 



/ Lion 

\4ivcE ion Un e 



1/2 



(2) 



= J^l(E)dE is the ionizing luminosity, and 
i/Qiom is the mean energy of the ionizing pho- 



where Li 

Eion — 

tons. 

Analyses of the broad emission lines in various AGN 
indicates that typical values of U for high ionization lines 
such as H/3 in AGN clouds are 0.1-1 and n e ~ 10 10 — 
10 11 cm~ 3 in the high excitation lines (cf. Rees, Netzer & 



Ferland 1989) so that Un - 10 
R ~ Q.037(n w UEi, 



9-11 



cm 



)-^ 2 L 1/2 vc 



(3) 

n e /10 10 cm- 3 and 



where L li5 = Lj OIi /10 45 erg/s, n w 
Ei on is in Rydbergs. 

This value (cf. Alexander & Netzer 1994), is in agree- 
ment with the results from recent reverberation mapping 
observations (e.g. Clavel et al. 1991, Peterson et al. 1991) 

c. Assuming the energy spectrum of the luminosity 
emerging from the central energy source can be approxi- 
mated by a power law 1(E) oc E~ a we find 



Eion, — 



1 



-En 



a > 1 



(4) 



(Since we are interested only in the ionizing spectrum, 
between Eq —1 Rydberg and a few keV, the power law 
assumption may be a reasonable approximation. Equation 
4 holds only for a > 1, as otherwise Li on and E diverge. 
This divergence is not physical, as actually the hard X-ray 
spectrum cuts off at at 50-100 keV). 

Combining Eqs. (1), (3) and (4) we have 



Fig. 2. X-ray continuum slope, obtained from a power-law 
fit, versus FWHM of H/3 for Seyfert 1 galaxies (taken from 
Walter & Fink 1993) and NLS1 (BBF96).The size of the circles 
corresponds to the ROSAT luminosity. 



FWHM « (UOOkm/s) (^^) ^ L^M^ 2 (5) 

where M 7 = M/10 7 M Q . Since we do not know the ioniz- 
ing luminosity, we try to express it in terms of the X-ray 
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luminosity and the spectral slope. The observed X-ray lu- 
minosity L x is in the ROSAT 0.1-2.4 kcV band, while the 
ionizing luminosity extends to lRyd=13.6eV. Assuming 
the power law spectrum in the ROSAT band may be ex- 
trapolated to a lower energy E we have for the ratio 



E {keV) 



l-a 



(O.ITEq) 



l-a 



(6), 



L x 0.1 1 -" - 2.4 1 -" 

where here and below E without units indicates E in 
Rydbergs. Substituting this in Eq. (5) we have 



FWHMk, (HOOfcm/s) 



/ ni Ua \ 
\a-lj 



1/4 



-1/4 



J x45 



(0.14S ) (a - 1)/4 M 7 1/2 



(7) 

where L x45 = L x /10 45 erg/s. 

d. We do not know the BH mass of individual objects, 
but we can use a mass-luminosity relation in a statistical 
sense, over a large enough sample of AGN, as we did with 
the ionization parameter above. We assume that the mass 
is roughly proportional to the luminosity, that is, that 
the Eddington ratio has a relatively narrow distribution, 
compared to the range of the luminosity distribution over 
the sample (e.g. Wandel and Yahil 1985). Defining , r\ x — 
L x /L Edd where L Edd = (1.3 I0 45 erg/ s)M 7 , we can then 
express the central mass in terms of the continuum X-ray 
luminosity, M 7 « L X 4 5 /r/ x . Substitutinging this relation 
in Eq. (6) we have 



FWHM w (U00km/s)r]- 



1/2 ( n w Ua 



1/4 



L 



1/4 
:i;45 



(0.14£ ) ( "~ 1)/4 



(8) 



Detemining E and rj x will give us a relation between 
the line width, the spectral index and the X-ray luminos- 
ity. 

The ratio L x /LE dd is not known, but may be esti- 
mated from theoretical and observational arguments. The 
bolomctic L/Lem is often estimated to be in the range 
0.1-1. Observational estimates of L op t/ LEdd and L x /L E dd 
are in the range 0.01-0.1 (e.g. Wandel and Yahil 1985, 
Wandel & Mushotzky 1986), which may justify a choice 
of r] x = 0.1. 

Combining the parameters in Eq. (8) as P = 
(n w U)^ 4 (Vx/0.1)^ 2 gives 



FWHM w (5600fcm/s)1.6" 



E 



a—l 



a 



1/4 



r 1/4 p 



(9) 



Note that choosing Eq = 1 is equivalent to extrapolat- 
ing the power law spectrum to lRyd, which may overes- 
timate the ionizing luminosity, (especially for large values 
of a) if the spectrum actually flattens at a higher energy. 
However, because of the weak dependence on E , this ef- 
fect is small: for example, if L x = 10 44 erg/s, p — 1 and 



a = 2 E = lRyd gives FWHM ~ 1500/cm/s, while for 
E Q = 3Ryd we get 1800 km/s. 

In Fig. 1 we overlay the analytic relation between the 
FWHM and a with the observed data from Fig. (8) of 
BBF96 for different values of L x . The general shape of the 
observed distribution and the analytic relation of Eq. (9) 
are in good agreement. (Note that the convergence of the 
theoretical curves at a photon index of 2 (corresponding 
to a — 1) for large FWHM is not physical; if a high en- 
ergy cutoff would have been introduced, the curves would 
extend also to lower values). 

3. Predictions and tests of the model 

3.1. Line width 

Using Eq. (8) or (9) we may predict the line width from 
the observed luminosity and the X-ray continuum slope. 
Figure 3a shows the predicted line width (model) versus 
the actually observed one FWHM (observed). Ideally, one 
would expect all the points to be on the diagonal. However, 
the scatter in the parameters r\ x and nil (which, in the 
analytic model we have assumed to have a single value 
for all objects) and eventualy in the spectral shape of the 
EUV, can easily explain the scatter in the plot; infact, it 
is surprising the scatter is not larger than a modest factor 
of 3. 

The ratio between the observed FWHM and the calcu- 
lated velocity is scattered about unity, and shows no sig- 
nificant dependence on the other observables (L x or a) , as 
is expected if the model gives the correct dependence on 
these observables. The ratio appears to be weakly anticor- 
related with the observed FWHM (Fig. 3b), in the sense 
that very narrow line objects give larger values for the pre- 
dicted/observed FWHM ratio, and vice verca. However, 
the correlation coeffient is p = —0.2 which is insignificant, 
certainly compared with the p = 0.7 correlation coefficient 
of the predicted vs. observed FWHM. This indicates that 
the model probably takes the main effects into account, 
but there is a systematic residual dependence on the line 
width. 

3.2. Luminosity- dependent L x /M 

What does this residual dependence tell us? If the basic 
assumptions made above are correct, the bias has to be 
sought in the parameters. If, for example, L x /M depended 
on FWHM, so that it decreased with increasing FWHM, 
this would reduce the effect, since for objects with larger 
FWHM Eq. (8) would give larger calculated values and 
vice versa. One explanation for such a dependence is that 
objects with narrow lines tend to have a smaller L/M ra- 
tio, that is, narrow line objects tend to have relatively 
smaller mass black holes. 

We now present a different, explanation, by showing 
that such a dependence does indeed follow from the well 
known phenomenological relation between the optical-UV 
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Fig. 3. Line width prediction from the model. Left panel (3a): Predicted (from Eq. 9) versus observed FWHM of H/3 for Seyfert 
1 galaxies (taken from Walter & Fink 1993) and NLS1 (BBF96). Right panel (3b): Ratio of predicted and observed FWHM 
versus observed FWHM of H/3. 



luminosity (L ) and X-ray luminosity in AGN. In trans- 
forming the black-hole mass to luminosity we have as- 
sumed that the Eddington ratio (~ L/M) is approxi- 
mately the same for all objects. Since we use the X-ray 
luminosity, we should actually take into account the ratio 
L x /M, or, L x /L. It has been shown that in large samples 
of AGN L x oc L0.75±o.05 (- Kriss 19g8j Mushotzky and Wan- 
del 1989; note however that there is a considerable scatter 
around this relation). If so, assuming L Q /M ~ const, im- 
plies (on average) 

Vx ex L x /M cx L-°- 25 ~ L-°- 33 . (10) 

We can plug this into the theoretical relation (Eq. (8) 
or (9)) to get a phenomenological relation, which gives a 
modified dependence on L x : 

FWHM oc £°' 2 V 1/2 oc L° X A . (11) 

In order to see whether this explains the residual de- 
pendence on FWHM (Fig. 3b), we should express the 
FWHM calculated/observed ratio in terms of FWHM in- 
stead of L x . Expressing L x in terms of FWHM in Eq. (8) 
and using Eq. (10) gives r] x oc FWHM~ S . Since the val- 
ues of FWHM calculated in Fig. 3b assume a fixed value 
for L x /M, the calculated/observed FWHM ratio will have 
a bias of 1/rj 

x FW H M ^'^ which agrees with the 

dependence seen in Fig. 3b. 
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Fig. 4. Luminosity prediction from the model. Upper limit of 
the predicted ROSAT luminosity vs. observed ROSAT lumi- 
nosity. We assume an Eddington ratio of r\ = 1 and a relation 
between ionizing and ROSAT luminosity of L ion — 7 a ~ 1 L x . 
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3.3. Luminosities 

Eq. (8) predicts that the X-ray luminosity of individual 
objects should increase in the direction normal to curves 
of constant L x in the FWHM-a plane (see Fig. 1). Fig. 
2 gives the actual observed ROSAT luminosities (the size 
of the symbols corresponds to L) indicating a general in- 
crease of the observed luminosities as predicted by the an- 
alytical relation. Reversing Eq. (9) and using the observed 
FWHM we have 



(pred) = 7° 



1 



En 



FWHM 
5600fcm/s 



(niot/rWO-l) 2 



(12) 



so we may calculate a predicted value for the X-ray 
(ROSAT band) luminosity, assuming e.g. rj x = 0.1 (Fig. 
4). Comparing this with the observed ROSAT luminos- 
ity may be an alternative test of the model. Theoreti- 
cally, if there were no scatter in the parameters, all ob- 
jects should lie on the diagonal. The scatter in Fig. 4 
appears larger than in Fig. 3a, since the parameters are 
in the fourth power. We also note that the slope of the 
best linear regression appears larger than unity, which re- 
flects the residual dependence seen in Fig. 3b. If we as- 
sume rj x does not depend on L x , Eqs. (11) and (12) give 
L x {pred) = L x (obs) 16 , approximately the right slope for 
the regression in Fig. 4. 

4. Determination of M from variability analyses 

In Sect. 2 (Eq. 8) we have assumed all objects have the 
same mass-luminosity ratio to derive an estimate for the 
black hole mass. In order to test this assumption we make 
below an alternative derivation with an independent esti- 
mate of the mass. 

An upper limit for the black hole mass is given by (e.g. 
Wandel & Mushotzky 1986) 



M(At) < (c 3 /10G)At ~ io 4 Ai M Q , 



(13) 



where At is in seconds. This relation assumes that the bulk 
of the X-ray continuum is emitted within 5 Schwarzschild 
radii, and that the variability is not affected by beaming or 
relativistic motions (cf. Boiler et al. 1997 for an example 
of relativistic boosting effects in a radio-quiet, ultrasoft 
NLS1). 

Using the doubling time as the characteristic time 
for variability we determine upper limits for the black 
hole mass for those objects of BBF96 and Wandel & 
Mushotzky (1986) whose continuum luminosities have 
been observed to vary significantly. Fig. 5 gives the dou- 
bling time (we have extrapolated amplitude variations lin- 
earely to a factor of 2 to determine At) versus the observed 
FWHM. Although there is some scatter, there is a strong 
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Fig. 5. Parameters determined from variability analyses. Dou- 
bling time versus FWHM of H/3 for NLS1 (filled circles) and 
normal Seyfert Is (open circles) . 



indication for an increase of the doubling time with ob- 
served FWHM. This may indicate that NLS1 tend to have 
lower-mass black holes than ordinary Seyfert 1 galaxies. 

Finally, for the objects with an established doubling 
time and mass upper limit, we may calculate the model ve- 
locity without having to assume any L/M relation. Com- 
bining Eqs. (7) and (12) we have 



FWHM < (3500fcm/s) 



At 
10 4 s 



1/2 



a - 



1/4 



L 



-1/4 



x45 



(14) 

Fig. 6a shows the predicted FWHM upper limits, and 
6b - the ratio FWHM(observed)/FWHM(predicted), vs. 
FWHM (observed) for the objects with doubling times. We 
see that the predicted FWHM is larger than the observed 
one by an average factor of 10, which is not surprising, 
considering that the variability method gives only an up- 
per limit for the mass. Similarly to the photoionization 
method, also here the ratio shows a slight anticorrelation 
with the observed FWHM. 



5. Summary and Discussion 

5.1. Basic Model 

We derive a simple theoretical model which explains the 
observed distribution of line width and spectral index for 
Seyfert 1 galaxies. The observed velocity dispersion is de- 
termined by the mass of the black hole and the broad 
line region radius. Sources with steep X-ray spectra (for a 
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Fig. 6. FWHM predicted from the model using X-ray variability analyses to set an upper limit on the black hole mass (Eq. (9)). 
Left Panel (6a): Prediced vs. observed FWHM. Right pannel (6b): FWHM(observed)/FWHM(predicted) vs. FWHM(observed) 
for the objects with doubling times (Eq. 14) 



given luminosity) have a stronger ionizing power than flat- 
spectrum sources, which implies that the BLR is formed 
at relatively larger distances from the central source, and 
hence has a smaller velocity dispersion and a smaller ob- 
served FWHM. 



case for the parameter r\ x — L x /LEdd, and is treated 
in Sect. 3.2. In principle, also the other parameters may 
systematically depend on some variable, but the data re- 
quired to check this (e.g. separate determinations of Un e ) 
for individual objects are not yet available. 



5.2. Model Parameters 

The model predicts a particular distribution of the lumi- 
nosities in the three dimensional FWHM-Lj — a space, 
which is essentially confirmed by the data. The FWHM 
also depends on the parameters U, n, Eq (and perhaps 
other parameters we are not aware of), as shown in Eqs. 
(8) and (9). All those parameters, however, cannot be ob- 
served directly or are not available for the individual ob- 
jects, so we develop and test the explicit teoretical relation 
for the observables that are available for the individual ob- 
jects in our sample, and show that the data is consistent 
with the theoretical relation we derive. The dependence 
on the other parameters is kept in the equations, but be- 
cause of our observational ignorance we assume that all 
objects have the same values for those parameters. This 
choise leads to an implicit error, which contributes to the 
scatter in Fig. 3. Because of the weak dependence on the 
parameters, the expected error is relatively small, of the 
order of 3 for the expected parameter range. 

In addition to the scatter in the parameters, there 
may be also a systematic dependence on some of the ob- 
servables, which would lead to a dependence of the pre- 
dicted/observed FWHM ratio. This is shown to be the 



5.3. Observational Bias 

BBF96 searched the literature for reports of NLS1 and 
found 46 in total. 32 of these NLS1 were located in in the 
fields of view of ROSAT pointings in the public archive. 
Only one of these 32 objects was not found above a 5 a 
detection limit. Walter & Fink 1993 selected all Seyfert 
1 type AGN with more than 300 counts in the ROSAT 
All-Sky Survey data which were observed at least once in 
the ultraviolet with IUE. Although the selection of broad 
and narrow-lined AGN is based on these special selection 
criteria, there are no obviousg selection effects that can 
account for the absence of broad line AGN with steep soft 
X-ray spectra. The measurements in the photon index and 
the Hg values are totally independent so that there is no 
obvious way a spurious correlation could be induced. 

5.4- Independent determination of the Mass 

Estimating an upper limit for the central mass directly 
from X-ray variability data we find that AGN with narrow 
optical emission lines may have lower black hole masses 
than broad line AGNs (Fig. 5). In a separate work we will 
consider physical parameters leading to soft X-ray spectra 
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to further improve our model. In order to test this mass 
estimate, it is possible to compare the BLR radius cal- 
culated from the ionizing luminosity (e.g. Eq. 2) or from 
reverberation mapping to the radius inferred from the ob- 
served FWHM and the mass estimated from the variabil- 
ity. Doing this indicates that the variability mass tends 
to be larger than the real mass for low narrow-line AGN 
(Wandel and Boiler 1997). 
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5.5. What causes steep X-ray spectra in NLSls? 

We have shown that sources with a steep soft X-ray spec- 
trum will have narrower broad lines than objects with flat 
spectra. A question we have not touched in this work is 
the cause for the observed steep X-ray spectra in many 
NLSls. In principle, smaller central masses could produce 
steep soft X-ray spectra as follows. For a given accretion 
rate, lower mass black holes could yield a hotter accretion 
disk, which would radiate more energy in the soft X-ray 
band. X-ray reprocessing by the accretion disk would pro- 
duce a soft X-ray spectrum (Matsuoka et al. 1990, Pounds 
et al. 1990) that extends to higher energies. This would 
make them appear to have steeper ROSAT spectra since 
there would be more flux in the soft X-ray band. This 
scenario will be considered in a separate work. 
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